General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



X-621 •75-254 

rr.crrviNt 


L- . 70^^^ 

NEUTRAL COMPOSITION 
GRAVITY WAVES IN THE THERMOSPHERE 
OBSERVED BY ESRO 4 


V“ W* V^’S T*I 

ny rc~0 


vrti >^pn. ^T"^T0N 

TH*: Tiicoxo-^rjH*'' " OQfiEFVTf' 
?* V Hr n. rnn " 




V7i>- 1 


linclan 

S 


H. TRINKS 
H. G. MAYR 


OCTOBER 1875 


80DDARD SPACE FLIGHT CENTER 

QREENRILT, MARYUND 


X-<i21.75-2r»4 


NKl'TRAI. ^OMP()SITU)N GILWITY WAVES IN TUI. EHMOSPHERE 

OBSERVED BY ESRO I 


II. TRINKS* 

II. G. MAYR 

Laboratory’ for Planetary Atmosphere 


♦NAS/NRC Resident Research Associate, on leave from 
Physilcdisches Institut, Universitaet Bonn, Germany 


GODDARD SPACE FLIGHT CENTER 
GREENBELT, MARYIJ^ND 


NllUTRAL COMPOSITION GRAVITY WAVT'.S IN THi; TIIIRNWSPHIRI-. 


OBSl'RVF.I) BY r.SRO 4 
ABSTRACT 

Neutral composition waves with a wave length of about 5000 km and 
a wave per'od of about 2.5 hours were observed by fSRO 4 in the altitude 
region of 250 km. The amplitudes arc of the order of 25*» for Ar, 15* 
for N^, and roughly 10*. for lie and 0 at 25" geographic latitude. The 
Ar and N\ waves are almost in phase, whereas He is in anti -phase and 0 
is in between. The wave amplitudes arc seen to decrease towarils lower 
latitude, uggesting that the composition waves are launched by auroral 
sources. iltancous ground based ionosonde miasurcments of the F2 

layer critical frequency at mid and low latitudes s'.ow a wave period con- 
sistent with the satellite observations. From the relation between wave 
length and wave period these waves are identified as gravity waves. The 
theoretical investigation with a multi-component model shows that diffu- 
sion plays a major role in explaining these wave phenomena. ITie phase 
and •'mplitude relation between atmospheric constituents are sensitive to 
the altitude region in which energy is deposited thus suggesting that 
.loule heating or soft particle precipitation are the predominant energy 
sources . 
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NIUTRAL COMPOSITION OP.*.v:,r KAVl.S TMi: THPRMOSI IIERF. 

OBSI-RVMI) JY I.SRO 4 

INTRODUa iON 

Since Mines [I960] postulated internal atmospheric gravity waves 
only a few observations of wave- like structures of the noi:ti-al gas density 
have been published (e.g. Newton et al ., 19f>9; Forbes a n d Marcos , 197S|. 
Recently, measurements in the lower and middle thermosphere obtained by 
satellite borne neutral gas mass spectrometers were reported by I led in 
ct ai . , (1974), Proe l ss and von Zahn (1975| revealing the existence of 
wave structures in the neutral composition, the "wavelength" ranging 
between about 30 km and several 100 km. This paper presents neutral 
composition data obtained by the gas analyzer aboard KSRO 4 during four 
consecutive perigee passes showing large scale wave structures which are 
inter^ireted as gravity waves. 

bXPLR lMPN TAl. MLTI lOD 

cSRO 1 was a spin- stabi 1 i zed polar orbiting spacecraft with an 
inclination of 91®. During the time interval to be discussed here the 
perigee height was 23.S km while the apogee height was .about 500 km, thus 
the height variation within the considered latitude range of 120“ was 
only 80 km. The perigee latitude was about 38“ S, the local solar time 
at perigee was 13:36. Tl»e gas analyzer [Trinks and von Z.ab n, 1975] 
measured the number densities of the total oxygen content n(0) ♦ 2 nfO^), 
molecular nitrogeii n(N.,), helium n(llc) and argon nfAr). The measured 
oxygen content reflects mainly the atomic oxygon density since the con- 
tribution of molecular oxygen becomes negligible under most circumstances 
above the F.SRO 4 perigee. 
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MITHOD OF DATA ;\NAI.YS1S 


For this study the density data of four orbits obtained between 
1:15 and 6:20 irr on (> March, 1974 within the geographic latitude range 
15*N. . . 89*S. . . 75®S were selected. From these data it was apparent that 
altitude, latitude and local solar time dependence of the density of each 
atmospheric constituent in most cases exceeded the wave induced struc- 
tures. Consciiuent ly , the variations not as.iOciated with wave structures 
had to be filtered out. 

To eliminate the altitude dependence the densities of all constituents 
except helium were converted into exospheric temperatures using Jacchia's 
model [.lacchia , 1971). The latitude and local solar time dependence of 
the obtained temferatures were filtered out by applying a least square 
analysis in terms of low order spherical harmonics. The local solar 
time dependence was included to distinguisli both sides of the south pole 
Siimpled by FSRO 4. In order to fit only those structures which re 
unchanged from one pass to the next, the fit for each constituent 
included the data of all four passes. Thus a temperatu.'e fit for 0, 
for N, and for Ar was obtained. The exospheric temperature (its were 
then reconverted into densities using again .Jacchia's model. In the case 
of helium the density ilata were directly fitted by spherical harmcnics 
since the altitude variations in the height range of 80 km were small. 
Figure 1 represents the resulting density fits versus geographic latitude; 
also shown is the latitude dependence of the satellite altitude during 
the passes considered. In the high.-r latitude region the neutral composi- 
tion obviously is disturbed (increased Ar and density, decreased He 


2 


and 0 doti5ity) in response to disturbed jjeomaijnet i c conditions; diirii.n 
this time interval a^^ changed from 2 ~ to 

In order to separate the wave structures the original density data 
were finally conmared with the density fits by computing the ratio 
R (W) for each constituent and orbit, P being the geographic latitude. 

Since the same fit was used for all passes, the centerlines of the 
R (P) data woie slightly different from unity which was especially true 
at higher latitudes. I’herefore the ratios R(W) were fitted by the 
spherical harmonics and Pj for each constituent and orbit. I'he old 
ratios were then divided by the corresponding fit values and replaced 
by thesi new ratios. This procedure elimi""ted the global sc.ile variations 
from one pass to the other without significantly effecting the wave 
structures. The results are presented in Tigure 2 for \r, N,, 0 and Me. 

In order to emphasize the wave structures with la “go wavelengths (which 
are to be investigated in this paper), smoothed curves were hand fitted 
to the data points neglecting small scale variations. 

DISaiS.S ION 

At low and midlatitudes a well pronounced sinusoidal wave stri 
develops in the case of Ar and N, ind weaker but still recognizable 
in the case of 0 and Me. The "wavelength" is of the order of .'>0“ lati- 
tude corresponding to about .SOno km which has to be considered as an upper 
limit of the real wavelength since the satellite's trajectory was not 
necessarily perpendicular to the wave front. The phases of the N, wave 
show a slight phase difference when compared with the Ar wave, whereas 
both 0 and Me are almost in antiphase. In the higher latitude region the 
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large scale wave structure is more difficult to identify because small 
scale disturbances and waves become predominant. There is a well pro- 
nounced localized disturbance between 60* and 70“ geographic latitude 
during pass *1 and somewhat weaker during pass which is of the kind 
described by Wu If -Mat hies ct al ., (1974| and might define the source 
region of the observed waves. Small scale wave structures are esp ‘daily 
developed during pass *4; note that for this case theie is only a slight 
phase difference between the 3 constituents Ar, and 0, (a corresponding 
wave structure in He cannot be recognized due to the scatter of data 
points ) . 

To demonstrate the time development of the wave structure tb? 
smoothed curves of all four passes (see Figure 2) are put into one frame 
in the case of argon (Figure 3). Although it cannot be excluded that 
there is a longitudinal variation of the wave structure we assume that 
most of the observed changes from one pass to the next reflect a time 
development; in the follow»;ig, this assumption wilt bo strongly supported. 

Since we have four consecutive passes showing wave structures with 
almost the same amplitudes it is possible to fit the 4 data points 
(see Figure 2) obtained at a fixed latitude during all passes by a 
cosine function versus time t; 

R (*4p,t) = a cos (u't ♦ a), t = t (0^) for passes 1-4 

The fit in turn yields the ;imilitude a, the frequency w, and the phase a 
for e.'ch constituent. In order to reduce the influence from the scatter 
of the data points especially in the case of He and Ar mean values 
R (P^,t) were used including 5 data points ^A0,t) closest to the 


4 


sclcctwJ latitude P . To linearize the fit problem ui was set constant 
when applying the least mean square fit method. Varying u> by steps the 
minimum of the standard deviation finally determined the best fit for w. 

In the case of Ar (and N,) a well pronounced minimum of the standard 
deviation was found for a wave period 2 f/w ■ 140... ISO min. For 0 and 
lie in most cases no clear minimum of the stamlard deviation could be 
identified because of the small amplitudes of tne 0 and lie wave and the 
scatter of data points (helium). In T.ahle I the obtained wave narameters 
are listed for a wave period of 110 and l.SO min at different geographic 
latitudes. At 25"S the Ar w.ive shows an amplitude of about 25*» which is 
significantly dei 'easing towards the equator; furthermore a phase differ- 
ence of II is found between 25“S and the equator, yielding again a "wave- 
length" of about 59® in latitude. For the results are almost the 
same concerning the phase, where.is the amt»litudc is less by about 10'p. 

I'or 0 and lie the results from t'able I show a nhase difference of roughly 
().7n in the case of 0 and n in the case of Me with respect to the Ar wave. 
The .unplitudes are relatively small and become comparable to the standard 
deviations. (No results are presented at 0® and 12.5®.S geographic lati- 
tude where the stmdard deviations exceed the amplitudes). 

The waves observed by the liSRO 4 gas analyzer are in good agreement 
with ionospheric waves as deduced from Fr layer critical frequency data. 
Figure 1 represents f^F2 from three Australian and one New Caiincan 
station located at geographic longitudes between FSRO 4 passes >*2 and ^ 
(see Figure 2). Although the time resolution is poor and the fit curves 
in some cases may be questionable, the overall data clearly e.xhibit 


5 


ionospheric waves with a period of about 2.5 hours. The wavelength turns 
out to he about 50* in latitude, considering that the ionospheric waves 
above Brisbane (28®S) and Vanimo (5*S) are almost in antiphase as marked 
by the dashed line at 15:10 liNfT (5:10 UT). The wave maxima, which are 
presumably associated, are identified by arrows, indicating the progression 
of a wave towards lower latitudes. 

A model for the excitation of composition gravity waves has recently 
been discussed by Mayr and V o Hand (!y75) where it was shown that wind 
induced diffusion is sigriflcant In producing phase differences between 
atmospheric constituents. Thi; equations of mass, energy and momentum 
conservation were solved for a multi -component ges (including N^, 0 and 
He) considering norizontai and vertical diffusion. Time and latitude 
dependence were separated in terms of Fourier and spherical harmonics; 
longitudinal and local time variations were ignored. The heat source 
responsible for exciting these waves was assumed to peak near 140 ^m due 
to .Joule heating in the auroral zones. The upper and lower boundaries 
of the model were at 700 and 80 km respectively, presumably removed far 
enough from the excitation region so that homogeneous boundary conditions 
are justified. For our interjiretat ion we adopted this model but extended 
it to include argon 

Adopting the spherical harmonic corresponding to an average wave 
length of 5000 km (in accordance with the predominant wave structure in 
our data) the thermcspheric response was evaluated for a scries of 
frequencies. ITie resulting ;implitudes and phases of Tg, 0 and S '2 are 
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shoton for 250 km in I'lgurc 5. It rovcals resonance maxima near 6 x 10 
see'* (t « 2.0 hoiirsl indicating that most of the v*ave energy for the 
5000 km w;!,c length is cairied by a wave period of ,\.0 hours, in reason- 
able agreement with our observations. ince the physical mech.anisms for 
prndui ing the resonance maximum also led to the dispersion relation in 
gravity wave theory we identify these waves as gravit> waves. The rapid 
phase transitions near the "resonatice frequency" are responsible for the 
wave propagation as was shown by Mayo and Vo 11 an d 11975'. By comparison, 
the phase differences between 0 and (and similarly between other 
constituents) remain relatively constant which makes then particularly 
useful in comparative studies regarding the excitation source for example. 

Tor I* and w ■ 6 x 10* * sec*^ the .amplitude and phase distributions 

O 

of tempe’.ature and composition arc shown in I'igure h. ITtc relative 
■amplitudes of the atmospheric constituents are in substanti.il agreement 
with the KSRO 4 data at 250 km. The extent to which deviations from 
diffusive equilibrium are important for N,, 0 and lie has been discussed 
by Mayr and Vo 1 land 11975). For \t the .amplitude of the "effective 
scale height temperature" (defined in May r e t a 1 . 119741) exceeds the gas 
temperature by 70*o at 150 km and .t0% at 250 kn, which contributes sig- 
nificantly to the comparatively large Ar density .amplitude. It is 
apparent from Figure 6 that the phase differences between individuil 
constituents agree qualitatively with the F.SRO 4 measurements at 250 km 
altitude: Ilic difference is small between and Ar ( 20*) and largest 

between Ar and He ( 110®); the phase ditference between Ar and 0 (60®) 

is in between. 
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Apart from the uncertainties in the phase determination of satellite 
data there are considerable uncertainties in the excitation of composi- 
tion waves. Momentum sources associated with ion convection, and local 
time or lont(itudinal variations could be important. Furthermore, the 
height structure of the heat source is not well determined, depending on 
the heighc distribution of the electron density in the case of Joule 
heating and depending on the energies of precipitating particles. 

Moving the heat input maximum from 140 km (the case discuised 
earlier) up to 200 km has pronounced effects on the phase relations 
between atmospher c const ituents. Iliis is illustrated in Figun 7 where 
the phase differences with respect to Ar arc shown as a function of 
height. Kith the exception of N, the Ar-0 and Ar-He phase differences 
arc considerably larger for the 200 kiu than for the 140 km heat source, 
jn better agreement with the FSRO I data. From the comparison between 
the experimental results and the theory we are led to conclude that 
thermospheric sources, possibly as high in altitude as 200 km, are 
responsible for the cxiitation of the observed wave phenomena. 
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Tabic 1. Period 2n/w, amplitude a, and phase o as derived from best 
fit R(W^,t) • a cos (u»t ♦!) at fixed KcoKraphlc latitudes R(0^,t) 
us in I'ig. 2. Also listed is the minimum standard deviation S.Dcv. and 
the halfwidth A(2s/w) ■ f2w/wj) - (2s/ui,), Wj and correspond to twice 
the minimum standard deviation (see also text). 

0 (las 2n/u) a <x S.OnV. A(2"/w) 

o 


(•S) (min) (\| (w) [\] (min) 


25.0 

Ar 

140 

25.5 

1.28 

2.8 

15 



150 

24.4 

1.11 

4.4 



140 

14.4 

1.27 

5.2 

30 



>50 

15.0 

1.12 

3.6 


0 

i -.0 

6.3 

2.0 

5. 1 




150 

7.3 

1.84 

5.0 



He 

140 

7.6 

2. If 

6 . 5 




150 

9.0 

2.06 

6.0 


12.5 

Ar 

140 

9.0 

1. 72 

1.3 

15 



150 

9. 1 

1.64 

0.6 



140 

4.5 

1.56 

5.7 




150 

4.5 

1.42 

5.4 


0 

Ar 

140 

8.8 

2.32 

'y 'y 

35 



150 

7.7 

2.15 

3.3 
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DENSITY (cm 3) ALTITUDE 



Figure 1 : Fits of atomic oxygen, molecular nitrogen, helium ami argon 

number densities obtained during 4 orbits (passes 1-4) between 
1:15 and 6:20 UT on 6 March, 1074. The uppermost curve rep- 
resents the altitude of observation. The local solar time at 
perigee was 1.5:56. The density fits are used to m; ke wave 
structures recognizable (see Fig. 2). 
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Figure 2 : iJcnsity ratios R equal to the measured nunber density data 

divided by the corresponding density fit value (Figure 1) for 
Ar, N*, 0 and He during four consecutive satellite passes. 

In order to eliminate variations of global scale occurring 
during the four considered passes the density fit for each 
constituent was allowed to vary by *’iC- spherical harmonics 
P and P before computing R for each pass. The ordinate is 
a logarithmic scale. The latitude region 90® *1® was not 
sampled by the F.SRO 4 (inclination 91®). The ^ hourly index 
a was decreasing from 22 to 7 for passes 1 to 4. The single 
data points were hand fitted by smoothed curves to emphasize 
the large scale wave structures. 
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Fig ure 3 ; Smoothed density ratios R for argon during four consecutive 
satellite passes from Figure 2. Notice the time development 
of the wave structure. 
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STATION 
GEOGR. POS. 

CANBERRA 
35 S. 149E 


BRISBANE 
28S, 153E 


TOWNSVILLE 
19S, 147E 


VANIMO 
3S, 141E 


Figure 4 . Time variation of the F2 layer critical fiequency of four 

ionosonde stations located at approximately the same geographic 
longitude which is between the longitudes of liSRO 4 passes 2 
and 7 > (see Fig. 2). The bar indicates the magnitude of a 
2 Mllz change in f F2. The wave maximum marked by an arrow 
at each station i§ believed to be the same maximum migrating 
from midlatitudes towards the equator. The triangles on top 
of the figure indicate the times of the four tSRO 4 passes 
at the geographic latitude of (UT»FNfr-10 hours). Note, 

that during the time indicated by the broken line the waves 
above Brisbane and Vanimo are in antiphase, correspending to 
a wavelength of about 5000 km. 
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Figure 5 : Computed relative amplitudes and phases of N- and 0 for 1’^^ 
as a function of frequency u at 250 km. Also shown is the 
phase di ffercnco jbetw^en 0 and N' . A resonance frequency occurs 
near ui » 6 x 10“’sec" corresponaing to a wave period of 2.9 
hours . 
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Figure 6: Computed relative amplitudes and phases of Ar, N_, 0, He and 
Tg as functions of height for P and u ■ 6 x TO'^sec"^. 
Note that the temperature is out of phase with N. and almost 
in phase with He. ^ 
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PHASE DIFFERENCES WITH RESPECT TO Ar (degrees) 


Figure IMiasc differences for N,. 0 and He with rcRcrd to \r computed 

~ for two different heat Sources, one peaking near HO km (cor- 

responding to the results in Figures 5 and (■>) and another one 
which peaks near 200 km. 
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